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ABSTRACT 

We present a study of optical Fe II emission in 302 AGNs selected from the SDSS. We group 
the strongest Fe II multiplets into three groups according to the lower term of the transition 
(b 4 F, a. 6 S and a 4 G terms). These correspond approximately to the blue, central, and red 
part respectively of the "iron shelf" around H/3. We calculate an Fe II template which takes 
into account transitions into these three terms and an additional group of lines, based on a 
reconstruction of the spectrum of I Zw 1. This Fe II template gives a more precise fit of the 
Fe II lines in broad-line AGNs than other templates. We extract Fe II, Ha, H/3, [O III] and [N II] 
emission parameters and investigate correlations between them. We find that Fe II lines probably 
originate in an Intermediate Line Region. We notice that the blue, red, and central parts of the 
iron shelf have different relative intensities in different objects. Their ratios depend on continuum 
luminosity, FWHM H/3, the velocity shift of Fe II, and the Ha/H/3 flux ratio. We examine the 
dependence of the well-known anti-correlation between the equivalent widths of Fe II and [O III] 
on continuum luminosity. We find that there is a Baldwin effect for [O III] but an inverse Baldwin 
effect for the Fe II emission. The [O III]/Fe II ratio thus decreases with Lasioo- Since the ratio is 
a major component of the Boroson and Green eigenvector 1, this implies a connection between 
the Baldwin effect and eigenvector 1, and could be connected with AGN evolution. We find that 
spectra are different for H/3 FWHMs greater and less than ^3000 kms -1 , and that there are 
different correlation coefficients between the parameters. 

Subject headings: atomic processes - galaxies: active - (galaxies:) quasars: emission lines 



1. Introduction 

Optical Fe II (AA 4400-5400 A) emission is one 
of the most interesting features in AGN spectra. 
The emission arises from numerous transitions of 
the complex Fe II ion. Fe II emission is seen in 
almost all type-1 AGN spectra and it is especially 
strong in narrow- line Seyfert Is (NLSls). It can 
also appear in the polarized flux of type-2 AGNs 
when a hidden BLR can be seen. Origin of the op- 
tical Fe II lines, the mechanisms of their excitation 
and location of the Fe II emission region in AGN, 
are still open questions. There are also many cor- 
relations between Fe II emission and other AGN 
properties which need a physical explanation. 



Some of the problems connected with Fe II 
emission are: 

1. There have been suggestions that Fe II 
emission cannot be explai ned with standard 



photoionization models ( see Collin fe Jolyl [2000; 
Collin-Souffrin et al.ll980l and references therein) . 



To explain the Fe II emission, additional mech- 
anisms of excitation ha ve been propos ed: a) 



continuum fluorescence |Philhpsl 11978a , 19791). 
b) collisional exc itation (jjolvl Il98ll Il987l Il99lt 
Veron et al ]|2006l). c) self- fluorescence among Fe II 



transitions ( Netzer fe Wills! Il983l) d) fl uorescent 
excitation by the Lyq and Ly/3 lines (IPenston 



19881 IVerner et ail ll999) . 
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2. There are uncertainties about where the 
Fe II emission line region is located in an AGN. 
A few proposed solutions are that the Fe II lines 
arise: a) i n the same emission reg ion as the broad 
H/3 line (|Boroson fe Green 19921) . b) in the ac- 



With increasing continuum luminosity, the EW of 
the C IV li nes decreases, ( t he so-called "Bald - 



cretion disk near the central black hole, which 
produces the double-peaked broad Balmer emis- 
sion l ines ( Collin- Souffrin et al. 198C)t Zhang et al.l 
2006), c) in a region which can be heated by 
shocks or from an inflow, located between BLR 
and NLR - the so-called the Intermediate Line 



Region (Marziani & Sulentic 1993; Popovic et al 
20041 120071 120091 : iHu et al.l l2008al iKuehn eTal 



20081 ) and d) and in the shi elded, neutral, outer 
region of a flattened BLR (jGaskell et al.l 12007 ; 
Gaskel]|l2009h . 

3. It has long been established that the 
Fe II emission is correlated with the radio, X 
and IR continuum. Fe II lines are stronger 
in spectra of Radio Quiet ( RQ) AGNs, than 
in Radio Loud (RL) ob jects (jOsterbrockl Il977 ; 
Bergeron fe Kunth 1984) . But, if we consider RL 
AGNs, Fe II emission is stronger in Core Domi- 
nant ( CD) AGNs than in Lobe Domin ant (LD) 
AGNs (ISetti fe Woltierlll977t [jolvlll99lh . The rel- 
ative strength of the Fe II lines correlates with 
the soft X-ray s lope, but anti-corr e lates with X- 



ray luminosity ( Wilkes et aL| Il987t iBoller et al. . 
19961 : lLawrence et al.l Il997l ) Strong optical Fe II 



emission is usually associated with strong I R lu- 
minosity ( Low et al. 19891 Lipari et al. 1993), but 
i t is anti-correlate d with IR color index a(60,25) 
(|Wang et al.ll2006l ). 

4. Many correlations have been observed be- 
tween Fe II and other emission lines. iBoroson fe Green 
(1992) give a number of correlations between the 
equivalent width (EW) of Fe II and properties of 
the [O III] and H/3 lines. The most interesting 
are the anti-correlations of EW Fe II vs. EW 
[O IIIJ/EW H/3, EW Fe II vs. peak [O III], EW 
Fe II/EW H/3 vs. peak [O III], and the correlation 
of EW Fe II/EW H/3 vs. H/3 asymmetry. An anti- 
correlation between EW Fe II and FWHM of H/3 
has been also found (lGaskellHl985t I Zheng fe Keel 
19911 : IBoroson fe Greenl I1992 K These correla- 
tions dominate the first eigenv ector of the prin- 
cipal component analysis of IBoroson fe Green 
( 19921 ) (hereinafter EV1). Correlations between 
EV1 and some properties of C IV A1549 A line s 
have also been pointed out (jBaskin fe Laori r2004). 



win effect", iBaldwinl Il977l IWang et al.l Il995h . 



while the blueshift of the C IV li ne and the 
EW of Fe II increase. Recently, iDong et~al 



(2009) lave found the Baldwin effect to be re- 
lated to EV1. Fe II lines are usually strong 



sorption Lines (BALs) 


Harting & Baldwin! 1 19861: 


Boroson & Mevers! 1992 


)• 



Searching for a physical cause of these cor- 
relations may help in understanding of the ori- 
gin of Fe II emission. Some of physical prop- 
erties which may influence these correlations 



are: E ddington ratio (L/LEd d) (|Boroson fe Green 
19921 iBaskin fe Laorl |2004| ). bl ack hole mass 
(Mrh) and inclination angle dMilev fe Miller! 
19791 ; [Wills fe Browndl 19861 iMarziani et aLlbOoli 

In addition to these properties, Fe II emis- 
sion strength may vary with evo l ution of AGNs. 
For example, Lipari fe Terlevich (2006!) have ex- 
plained some properties of AGN (as BALs, Fe II 
strength, radio properties, emission line width, lu- 
minosity of narrow lines, etc.) by an "evolutive 
unification model" . In this model, accretion arises 
from the interaction between nuclear starbursts 
and the supermassive black hole. This would mean 
that not only the orientation, but also the evo- 
lutionary state of an AGN, influence its spectral 
properties. Thus, young AGN have strong Fe II, 
BALs, weak radio emission, the NLR is compact 
and faint, and broad lines are relatively narrow. 
In contrast with this, old AGNs have weak Fe II, 
no BALs, strong radio emission with extended ra- 
dio lobes, the NLR is extended and bright, and 
the broad lines have greater velocity widths. 

In this paper, we investigate the Fe II emitting 
region by analyzing the correlations between the 
optical Fe II lines and the other emission lines in 
a sample of 302 AGN from the SDSS. To do this 
we construct an Fe II template. The strongest 
Fe II multiplets within the AA4400-5500 A range 
are sorted into three groups, according to the 
lower terms of the transitions: 3d 6 ( 3 F2)4s A F, 
3d 5 4s 2 6 5 and 3d 6 ( 3 G)4s 4 G (which we will here- 
inafter refer to as the F, S and G groups of lines). 
The F group mainly contains the lines from the 
Fe II multiplets 37 and 38 and describes the blue 
part of the Fe II shelf relative to H/3. The S group 
of lines describes the part of Fe II emission un- 
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der H/3 and [O III] and contains lines from mul- 
tiples 41, 42 and 43 Fe II. Finally, the G group 
contains lines from multiplets 48 and 49, and de- 
scribes the red part of the Fe II shelf. A simpli- 
fied Grotrian diagram of these transitions is shown 
in Figure [T] We analyze separately their relation- 
ships with other lines in AGN spectra (H/3, [O III], 
Ha, [N II]). In this way, we try to connect de- 
tails of the Fe II emission with physical properties 
of the AGN emission regions in which these lines 
arise. Each of the Fe II multiplet groups has spe- 
cific characteristics, which may be reflected in a 
different percent of correlations with other lines, 
and give us more information about the complex 
Fe II emission region. 

The paper is organized as follows: In Sect. 2 we 
describe the procedure of sample selection and de- 
tails of analysis; our results are presented in Sect. 
3; in Sect. 4 we give discussion of obtained results, 
and finally in Sect. 5 we outline our conclusions. 

2. The sample and analysis 
2.1. The AGN sample 

Spectra for our d ata sample are taken fr om the 
the 7 th data release ( Abazajian et al. 20091 ) of the 
Sloan Digital Sky Survey (SDSS). For the purposes 
of the work, we chose AGNs with the following 
characteristics: 

1. high signal to noise ratio (S/N > 20), 

2. good pixel quality (profiles are not affected 
by distortions due to bad pixels on the sen- 
sors, the presence of strong foreground or 
background sources), 

3. high redshift confidence (zConf>0.95) and 
with z < 0.7 in order to cover the optical 
Fe II lines around H/3 and [O III] lines, 

4. negligible contribution from the stellar com- 
ponent. We controlled for this by having the 
EWs of typical absorption lines be less then 
1 A (EW CaK 3934 A, Mg 5177 A and US 
4102 A>-1)0 Because of this it was not nec- 
essary to remove a host galaxy starlight con- 
tribution. 



5. presence of the narrow [O III] and the broad 
H/3 component (FWHM H/3 > 1000 kms" 1 ). 

We found 497 AGNs using the mentioned cri- 
teria. However, on inspection of the spectra we 
found that in some of them the Fe II line emis- 
sion is within the level of noise and could not be 
properly fitted, so we excluded these objects from 
further analysis (see Appendix A). As result, our 
final sample contains 302 spectra, from which 137 
have all Balmer lines, and the rest are without 
the Ha line (due to cosmological redshifts). The 
sample distribution by continuum luminosity (A 
5100 A) and by cosmological redshift is given in 
Figure [2j Corrections for Galactic extinction were 
made using an empirical selective extinction func- 
tion computed for each spectrum on the basis of 
Galac tic extinction coefficients given by ISchlegel 
(jl998l) and available from the NASA/IPAC Ex- 
tragalactic Databas^E Finally, all spectra were 
de-redshifted. 

To subtract the continuum, we used the DIPSO 
software, finding the continuum level by us- 
ing continuum windows g iven in the paper of 



Kuraszkiewicz et al.l (|2002l ) (see Figure [3]). The 
used continuum windows are: 3010-3040 A, 3240- 
3270 A, 3790-3810 A, 4210-4230 A, 5080-5100 A, 
5600-5630 A and 5970-6000 A. We used the same 
procedure to subtract the continuum in all ob- 
jects. Continuum subtraction may cause systemic 
errors in the estimation of line parameters and 
continuum luminosity. We estimated that the er- 
ror bars are smaller than 5%, especially within the 
Ha+[N II] region and in the red part of the "iron 
shelf". 

We considered two spectral ranges: AA 4400- 
5500 A and AA 6400-6800 A. In the first range, 
dominant lines are the numerous Fe II lines, [O III] 
AA4959, 5007 A, H/3 and He II A4686 A, and in the 
second range, Ha and [N II] AA6548, 6583 A. 

To investigate correlations between the Fe II 
multiplets and other lines in spectra, we fit all 
considered lines with Gaussians. The optical Fe II 
lines were fitted with a template calculated as de- 
scribed in the next section. 



2 http : / / ned www .ipac.caltech.edu/ 



1 For some objects we were able to find estimates of the stel- 
lar component in the literature and to confirm that the host 
galaxy fraction is indeed small in our sample (see Appendix 
A). 
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2.2. The Fe II (AA 4400-5500 A) template 

A number of authors have created an Fe II 
template in the UV and optical range (see 
Vestergaard fe Wilkesl2001 , and references therein) . 
Boroson fe Greenf l 1992h applied empirical tem- 
plate by removing all lines which ar e not Fe II, 

from t he spectrum of I Zw 1. Similarly. I Veron-Cettv et"al 
(|2004[ ) constructed an Fe II template by identify- 
ing systems of broad and of narrow Fe II lines in 
the I Zw 1 s pectrum, and meas uring their relative 
intensities. iDong et al.l fioOS) improve on that 
template by using two parameters of intensity - 
one for the broad, one for the narrow Fe II lines. 
All these empirical templates are defined by the 
line width and the line strength, implying that rel- 
ative strengths of the lines in the Fe II multiplets 
are the same in all objects. 

In theoretical modeling, significant effort has 
been made to calculate the iron emission by 
including a large number of Fe II atom i c lev- 
els, going to high energy (fVerner et aT 



Sigut fe Pradhanll2003l:ICollin-Souffrin et al 



1999: 



1980: 



Bruh weiler fe Vernerl 120081 ). iBruhweiler fe Verner 
(|2008l) calculated Fe II emission using CLOUDY 
code and 830 level Fe II model. In this model, en- 
ergies go up to 14.6 eV, producing 344,035 atomic 
transitions. For their calculations they used so- 
lar abundances for a range of physical conditions 
such as the flux of ionizing photons [$h] , hydrogen 
density [na], and microturbulence [£]. 

Using existing Fe II templates, we found that 
empirical and theoretical models can generally fit 
NLSyl Fe II lines well, but that in some cases of 
spectra with broader H/3 lines, the existing models 
did not provide as good fit. 

One of the problems in the analysis of Fe II 
emission is that it consists of numerous overlap- 
ping lines. This makes the identification and de- 
termination of relative intensities very difficult. 
Therefore, the list of Fe II lines used for the 
fit of the Fe II emission, as well as their rel- 
ative intensities, are different in different mod- 



els! Yjeron 1 Cetty_eLaL 20C 4: IBruhweiler fe Vernerl 



2008). Also, there is significant disagreement in 
values of osci llator strengths in different atomic 
data sources jFuhr et aL 1981 : Giridhar fe Ferrol 
19951 lKurucdll990h . 

To investigate the Fe II emission, we made an 
Fe II template taking into account following: (a) 



majority of multiplets dominant in the optical part 
(AA 4400-5500 A), whose lines can be clearly iden- 
tified in AGN spectra, have one of the three spe- 
cific lower terms of their transitions: F, S or G, 
and (b) beside these lines there are also lines whose 
origin is not well known but which presumably 
originate from higher levels. 

We constructed an Fe II template consisting of 
50 Fe II emission lines, identified as the strongest 
within AA 4400-5500 A range. 35 of them are 
sorted into three line groups according to the lower 
term of their transition (F, S and G). The F 
group consists of 19 lines (mainly multiplets 37 
and 38) and dominates in the blue shelf of the iron 
template (4400-4700 A). The S group consists of 5 
lines (multiplets 41, 42 and 43) and describes the 
Fe II emission covering the [O III] and H/3 region 
of the spectrum and some of emission from the 
red Fe II bump (5150-5400 A), and the G lines (11 
lines from multiplets 48 and 49) dominate in the 
red bump (5150-5400 A) (see Figure [Q). 

We assume that the profiles of each of lines can 
be represented by a Gaussian, described by width 
(W), shift (djU and intensity (I). Since all Fe II 
lines from the template probably originate in the 
same region, with the same kinematical properties, 
values of d and W are assumed to be the same for 
all Fe II lines in the case of one AGN. 

Since the population of the lower term is in- 
fluenced by transition probabilities and excitation 
temperature, we assumed that the relative inten- 
sities between lines within a line group (F, S or 
G) can be obtained approximately from: 



h. = il.il. e -{Ei-E 2 )/kT 

h Ai fa 52 



(1) 



where 1\ and I2 are the intensities of lines with the 
same lower term, Ai and A2 are the wavelengths of 
the transition, g\ and gi are the statistical weights 
for the upper energy levels, /1 and /2 are the os- 
cillator strengths, E\ and E2 are the energies of 
the upper levels of transitions, k is the Boltzman 
constant, and T is the excitation temperature. De- 
tails can be found in Appendix C. The excitation 
temperature is the sam e for all tran sitions in the 
case of partial LTE (see lGrieinlll997 ). but as it is 



3 Here we used W 



ay/2 



2^, d = 4^, where W D 
the Doppler width, AA is the shift and a is the velocity 
dispersion 
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shown in Appendix C, Eq (1) may be used also in 
some non-LTE cases. 

Lines from the three above mentioned groups 
explain about 75% of the Fe II emission in the 
observed range (4400-5500 A), but about 25% of 
Fe II emission cannot be explained with permitted 
lines for which the excitation energies are close to 
the lines of the three groups. The missing Fe II 
emission is around -4450 A, —4630 A, -5130 A 
and -5370 A. 

There are some indications that fluorescence 
proces ses may have a role i n producing some Fe II 



lines (jVerner et al.l I1999L Hartman fe Johansson 
2000). Process like self-fluorescence, continuum- 



fluorescence or Lya and Ly/3 pumping could sup- 
ply enough energy to excite the Fe II lines with 
high energy of excitation, which could be one of 
the explanation for emission within these wave- 
length regions. To complete the template for ap- 
proximately of 25% missing Fe II flux, we selected 
15 lines from the Kurutz database^ with wave- 
lengths close to those of the extra emission, up- 
per level excitation energies of up to —13 eV, and 
strong oscillator strengths. We measured their rel- 
ative inte nsities in I Zw 1 which ha s a well-studied 
spectrum ([Veron-Cettv et al 1 l2004h spectrum with 
strong and narrow Fe II lines. Relative intensities 
of these 15 lines were obtained by making the best 
fit of the I Zw 1 spectrum with the Fe II lines from 
the F, S, and G line groups. The extra lines are 
represented in Figure [1] (bottom), with a dotted 
line. 

Our template of Fe II is described by 7 free 
fitting parameters: width, shift, four parameters 
of intensity (for the F, S and G line groups and 
for the lines with relative intensities obtained from 
I Zw 1). The seventh parameter is the excitation 
temperature included in the calculation of relative 
intensities within F, S and G line groups. 

We found that for the majority of objects, the 
excitation temperature obtained from our fit is 
within the range: 9000 K - 11000 K (see Figure 
[p, which agrees we ll with theoretical predictions 
ICollin fe JohfeOOCh . However, as it could be seen 
in Eq. (1), our fit is not very sensitive to the 
temperature, especially for T>8000 K. Eq. (1) is 
also very approximate so estimated temperatures 



4 http: / / www.pmp.uni-hannover.de / cgi- 
bin/ssi/ test /kurucz / sckur . html 



should be treated with caution. 

To estimate the error in the excitation temper- 
ature, we found our best fit for a number of ob- 
jects. We then changed only the excitation tem- 
perature while keeping the other fit parameters 
fixed. We thus found Xmax f° r which the fit is 
still reasonable. After that, we constructed x 2 
vs. T and measured the dispersion in tempera- 
ture at 0-l-(Xmax~Xmin)- We obtained error-bars 
that were mostly within the range ±25%, with ten- 
dency for higher temperature to have a greater un- 
certainty. In further analysis we will therefore not 
consider the temperatures obtained from the fit. 

The wavelengths of the 50 template lines are 
presented in Tables Q] and [5J For the 35 lines from 
the three line groups we give multiplet names, 
configurations of the transitions generating those 
lines, oscillator strengths and the relative inten- 
sities. The intensities have been calculated using 
the Eq. (1) for excitation temperatures T = 5000, 
10000 and 15000 K (Table [Q. For the 15 lines for 
which relative intensities were measured in I Zw 1, 
we give wavelengths, configurations and oscillator 
strengths taken from Kurucz database, as well as 
their measured relative intensities (Tabled. 

The purpose of dividing the Fe II emission into 
four groups is to investigate the correlations be- 
tween the iron lines which arise from different mul- 
tiplets and some spectral properties, and to com- 
pare it with results obtained for the total Fe II 
within the observed range. Also, a template which 
has more parameters may give a better fit, since 
relative intensities among Fe II lines can be dif- 
ferent even in the similar AGN spectra, as is, for 
example, t he case for spe ctra of I Zw 1 and Mrk 
42 (as e.g. |Phillipsl[l978rA . 

We applied this template to our sample of 302 
AGNs from SDSS database, and we compare it 
with other theoretical and empirical models (see 
Appendix B). We found that the template fit Fe II 
emission very well (Figure [23] and Figure |2~41 . In 
the cases when the Fe II emission of an object has 
different properties from I Zw 1, small disagree- 
ments are noticed in the lines for which the rela- 
tive intensity was obtained from I Zw 1. However, 
the Fe II emission in these objects shows larger 
disagreement with the empirical and theoretical 
models considered in Appendix B (see Figure [25]) . 
Although our Fe II template has more free param- 
eters than other templates used for fitting the Fe II 
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emission, it enables more precise fit of Fe II lines in 
some AGNs (especially with very broad H/3 line) 
than the other Fe II templates we considered. 

2.3. Emission line decomposition and fit- 
ting procedure 

We assume that broad emis sion lines arise in 
two or more emission reg ions (Brothe rton et al 



1994 IPopovic et all 12004 iBon et al l 120061 12009; 



Ilic et al lteOolT so that their profiles are sums of 



Gaussians with different shifts, widths and inten- 
sities, which reflect the physical conditions of the 
emitting regions where the components arise. 

For this reason, we fit the emission lines we con- 
sidered ([O III] AA4959, 5007 A, R/3, He II A4686 
A, Ha and [N II] AA6548, 6583 A) with a sum of 
Gaussians, describing one Gaussian with 3 param- 
eters (width, intensity and shift from transition 
wavelength). 

Both lines of the [O III] AA4959, 5007 A dou- 
blet originate from the same lower energy level and 
both have a negligible optical depth since the tran- 
sitions are strongly forbidden. Taking this into ac- 
count, we assumed that the [O III] A4959 A and 
[O III] A5007 A lines have the same emission-line 
profile. We fit each line of the doublet with either 
a single Gaussian, or, in the case of a significant 
asymmetry, with two Gausssians. The Gaussian 
that describes the [O III] A4959 A, has the same 
width and shift, as the one which describes [O III] 
A5007 A line, and we took their intensity ratio to 



be 2.99 (|Dimitriievic et all 120071 ). [N II] AA6548, 



6583 A were fitted using the same procedure, with 
an intensity ratio of the doublet components of 
« 3. The He II A4686 A line was fitted with one 
broad Gaussian. 

To fit the Balmer lines a number of Gaussian 
functions were used for each line. It is usually as- 
sumed that there are two components for Ha and 
H/3: a narrow component representing the NLR, 
and a broad one representing the BLR. However, 
the broad emission lines (BELs), are very com- 
plex and cannot be properly explained by sin- 
gle Gaussian (which would indicate an isotropic 



(Brotherton et al.l 1994 


; Corbin & Boroson 19961 


PoDOvicetal. 2004: 


I 


ic et alJ 20061 Bon et al. 


2006L 120091: iHu et al. 




2008b) where the broad 



lines are assumed to be emitted from two kine- 



maticaly different regions: a "Very Broad Line 
Region" (VBLR), and an "Intermediate Line Re- 
gion" (ILR). We therefore tested fitting the broad 
lines with one and two Gaussians. We found that 
in most AGNs in our sample, the fit with two 
Gaussians was significantly better than the fit 
with one Gaussian. 

To apply a uniform model of the line fitting pro- 
cedure, we assumed that Balmer lines have three 
components from the NLR, ILR and VBLR. We 
exclude 4 spectra from the sample in which Ha 
line could not be decomposed in this way, so fi- 
nally our sample contains 133 spectra (from 137) 
with all Balmer lines which could be analyzed 
within the AA 6400-6800 A range. To check this 
assumption on the rest of the spectra, we exam- 
ined the corresponding correlations and we found 
that the luminosities and widths of the NLR, ILR 
and VBLR components of Ha are highly corre- 
lated with the same parameters of the NLR, ILR 
and VBLR components of H/3. This is in favor 
of our three component decomposition (Table [3]) . 
For the 133 objects, which have both H/3 and Ha 
lines, it can be noticed that 16 spectra have large 
redshifts of the VBLR H/3 component relative to 
VBLR of Ha, which causes disagreement and re- 
duces the correlation between the VBLR shifts 
of the rest of objects. Without these 16 objects, 
H/3 and Ha VBLR components correlate in width 
and shift, with coefficient of correlation ~ 0.40, 
P < 0.0001. An example of shifts between the 
Ha and H/3 VBLR components is shown in Figure 
03 The AGN SDSS J111603. 13 + 020852.2 has its 
VBLR H/3 component redshifted by «4258 kms -1 , 
while shift of the Ha VBLR component is «-1200 
knis . 

We assume that all narrow lines (and narrow 
components of the broad lines), originate from 
the same NLR, and thus expect that these lines 
will have the same shifts and widths. We there- 
fore took the Gaussian parameters of the shift and 
width of [O III], [N II], and the NLR components 
of H/3 and Ha to have the same values. 

We separately fit lines from the AA 4400-5500 
A range (Fe II lines, [O III], He II and H/3) with 
23 free parameters, and lines from the AA 6400- 
6800 A range (Ha and [N II]) with 8 free pa- 
rameters. We applie d a \ 2 minimization routine 
(jPopovic et al.ll2004j ). to obtain the best fit. Ex- 
amples of an AGN spectrum fit in both ranges are 
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shown in Figure [6] and Figure [7] 

2.4. Line parameters 

We compared the shifts, widths, equivalent 
widths, and luminosities of the lines. The shifts 
and widths were obtained directly from the fit (pa- 
rameters of width and shift). Luminosi ties were 
calcul ated using the formulae given in Peebles I 
(|1993l ). with adopted cosmological parameters: 
n M = 0.27, n A = 0.73 and tt k = 0. We adopt 
a Hubble constant, H Q = 71 kms _1 Mpc _1 . The 
continuum luminosity was obtained from the av- 
erage value of continuum flux measured between 
5100 A (AA5095-5100 A and AA5100-5105 A). 
Equivalent widths were measured with respect to 
the continuum below the lines. The continuum 
was estimated by subtracting all fitted lines from 
the observed spectra. Then, the line of interest 
(or the Fe II template), which is fit separately, 
is added to the continuum. Equivalent widths 
are obtained by normalizing the lines on contin- 
uum level, and measuring their fluxes (see Figure 
|8]) . To determine the full width at half maximum 
(FWHM) of the broad component of the Balmer 
lines, the VBLR and ILR components from the 
fit are considered as it is shown in Figure [9] The 
broad Balmer line (VBLR+ILR) is then normal- 
ized on unity, and the FWHM is measured at half 
of the maximum intensity. 

3. Results 

3.1. The ratios of Fe II line groups vs. 
other spectral properties 

Some objects have relative intensities of Fe II 
lines similar to I Zw 1, while others show sig- 
nificant disagreement, which usually shows up as 
stronger iron emission in the blue bump (mainly F 
lines) than one in the red bump (mainly G lines), 
comparing to I Zw 1. That objects could be fitted 
well with the template model which assumes that 
relative intensities between the Fe II lines from dif- 
ferent line groups (F, S and G) are different for 
different objects. 

Difference between the relative intensities of the 
Fe II lines from different parts of iron shelf (blue, 
red and central) can be illustrated well with his- 
tograms of distribution of the Fe II group ratio 
(F/G, F/S and G/S) for the sample (Figure [TOJ . 
The average value of the F/G luminosity ratio is 



1.44±0.55, of the F/S ratio 1.97±1.05, and of the 
G/S ratio 1.39±0.4€@. We found that the ma- 
jority of objects have Fe II group ratios close to 
the average values, but still significant number of 
objects are showing difference since the intensity 
ratios of Fe II line groups vary in different objects 
(see Figure [TO]) . The ratios of the Fe II line groups 
in I Zw 1 are indicated by the vertical dashed lines 
(F/G = 1.10, F/S = 1.56 and G/S = 1.41). 

We analyzed the flux ratios of the Fe II mul- 
tiplct groups and the total Fe II flux (see his- 
tograms in Figure lllj) . and we found that lines 
from multiplets 37 and 38 (group F) generally 
have the largest contribution to the total Fe II 
emission within the AA4400-5500 A range (about 
32%), while group S contributes about 18% and 
G about 23%. As expected, there are very strong 
correlations between the equivalent widths of Fe II 
line groups (Tabled]). 

Since line intensities and their ratios are indi- 
cators of physical properties of the plasma where 
those lines arise, we have investigated relations 
among the ratios of Fe II line groups with vari- 
ous spectral properties. A correlation between the 
FWHM H/3 and t he EW Fe II emission has been 
found previously (lGaskelll ll98"5 Boroson fe Greenl 



1992c ISulentic et al.ll2009t ). Consequently, one can 



expect that FWHM H/3 may be connected with 
the Fe II line group ratios. We therefore investi- 
gated correlations between the ratios of the Fe II 
groups (F/G, F/S, and G/S) and H/3 FWHM. 
We excluded the three outliers, which have neg- 
ligible flux for S group. As it can be seen in 
Figure [12] (first panel) there are different trends 
for spectra with FWHM H/3 less than and greater 
than ^3000 kms -1 . Therefore, we divided the 
sample into two sub-samples: 158 objects with 
FWHM B.j3 <3000 kms" 1 and 141 objects with 
FWHM H/3 >3000 kms" 1 (see Figure [HI and Ta- 
ble [5]). A simila r subdivision of AGN sam ples 
was performed bv lSulentic et al.l (|2002l . l2009l ). but 
for FWHM H/3 w 4000 kms -1 as the limiting 
FWHM of the subdivision. We found that FWHM 
H/3 » 3000 kms -1 is more a appropriate place 
to divide the objects into two groups, since all 
objects with strongly redshifted VBLR H/3 com- 
ponent (see Figure [5J belo ng to the FWHM H/3 



>3000 kms -1 subsample. ISulentic et ail (|2002h 



Here we give the averaged ratio value and dispersion. 
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also noticed that AGN with FWHM H/3 >4000 
kms -1 , have a larger redshifted VBLR H/3 (^5000 
kms -1 ). 

In Table [5] we present the correlations between 
the luminosity ratios of the Fe II line groups and 
other spectral properties such as: FWHM H/3, 
Doppler width of H/3 broad components, width 
and shift of iron lines and continuum luminosity 
(Lasioo)- They are examined for the total sam- 
ple and for two sub-samples, divided according to 
FWHM H/3, and denoted in Table H as (1) for 
FWHM H/3 <3000 kms" 1 and (2) for H^ FWHM 
>3000 kms- 1 . 

We found a difference in the correlations for 
those two sub-samples. For the H/3 FWHM > 3000 
kms -1 sub-sample (2) all three ratios (F/G, F/S 
and G/S) increase as H/3 width increases. On the 
other hand, for the H/3 FWHM < 3000 kms" 1 sub- 
sample no correlation is observed between these 
properties. 

The observed correlation between the G/S 
group ratio and H/3 FWHM may be caused by 
intrinsic reddening since reddening increase s as 
objects get more edge-on (IGaskell et al.ll2004h and 
the H /3 FWHM also increases (|Wills fe Browne! 
Il986l ). However, this cannot explain the F/G 
and F/S correlation with H/3 FWHM, and also, a 
strong influence of reddening on the G/S ratio is 
not expected because of the close wavelengths of 
lines from those groups. 

We also studied correlations between contin- 
uum luminosity (Lasioo) an d Fe II group ratios, 
since these may indicate excitation processes in 
the emitting region. We found that the F/G and 
F/S ratios decrease in objects where the contin- 
uum level is higher. Observed correlations are 
stronger for H/3 FWHM < 3000 kms -1 sub-sample 
(F/G vs. L A5 io : r = -0.51, P < 0.0001 and F/S 
vs. L A5 ioo: r = -0.41, P < 0.0001). The ratio 
of G/S seems not to be dependent on continuum 
luminosity. 

These correlations are the opposite of what 
would be expected from reddenin g, since re d- 



denin g decreases with luminosity (jGaskell et al 
2004). There are a few effects that could destroy 



expected reddening effect: (i) incorrect continuum 
subtraction, (ii) host galaxy fraction which de- 
pends on the Eddington ratio and (iii) different 
contribution of starlight in SDSS fiber which de- 



pends on the luminosity of the host galaxy and the 
redshift (see IGaskell and Kormendvl 2009h . There 
is a small possibility that these effects are impor- 
tant. As we already noted, the error-bar in the 
continuum subtraction is smaller than 5%. Also, 
we examined if the host galaxy contribution have 
influence on considered correlations. We found no 
correlations between the host galaxy fraction (de- 
termined in A5100 A) and Fe II group ratios@ The 
correlations appear instead to be connected with 
a weaker Baldwin effect for F group (see Table [7]). 

We also investigated the correlation of the log(L 
Ha/L H/3) with ratios of Fe II groups. We found 
correlations with F/G (r = -0.36, P < 0.0001) 
and with F/S ratio (r = -0.34, P < 0.0001), 
but there is no correlation with G/S (Figure fl~3l 
Table [6]), i.e., between ratio of red and central 
part. There is possibility that these correlations 
are caused by intrinsic reddening. 

3.2. Connection between kinematical prop- 
erties of Fe II lines and Balmer lines 
(Ha and H/3) 

We assume that broadening of the lines arise 
by Doppler effect caused by random velocities of 
emission clouds, and that shifts of lines are a con- 
sequence of systemic motions of the emitting gas. 
We therefore studied the kinematical connection 
among emission regions by analyzing relationships 
between their widths and shifts obtained from our 
best fits. 

The Doppler widths of Fe II and Balmer lines 
(H/3 and Ha) are compared in Figure [14J On X- 
axis we present Fe II width and on the Y-axis 
the widths of the H/3 (first panel) and Ha (second 
panel) components. The widths of the NLR com- 
ponents are denoted by triangles, the ILR com- 
ponents with circles, and the VBLR components 
with squares. Dotted lines show the average val- 
ues of the widths: vertical lines for Fe II compo- 
nents and horizontal lines for H/3 (or Ha). For 
the sample of 302 AGN (first panel) the average 
value of Fe II width and dispersion of the sample 
are 1430 ± 440 kms , while the average values 
for H/3 components are: 300 ± 150 kms -1 (NLR), 
1570 ± 700 kms" 1 (ILR) and 4360 ± 1440 kms" 1 



6 The amount of stellar continuum is taken from paper 
IVanden Berk et al.1 ll2006h for 106 common objects (see Ap- 
pendix A). 
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(VBLR). For the selected sample of 133 AGNs 
which have the Ha line in the spectra (second 
panel), the averaged value of the Fe II width is 
1310 ± 410 kms -1 and for the Ha components: 
240 ± 90 kms" 1 (NLR), 1160 ± 560 kms" 1 (ILR) 
and 4060±1650 kms" 1 (VBLR). It is obvious that 
the averaged Fe II width (~ 1400 kms -1 ) is very 
close to the average widths of the ILR H/3 and 
Ha components (1568 kms^ 1 and 1156 kms -1 re- 
spectively), while the averaged widths of NLR and 
VBLR components are significantly different from 
the Fe II width. 

Relationships among the widths of the Fe II 
and ILR components are also presented in Figure 
[TBI (Table EJ). The correlation between the Fe II 
width and the width of Ha ILR is r = 0.77, P < 
0.0001, and between the Fe II width and the H/3 
ILR width it is r = 0.73, P < 0.0001. We also 
found a correlation between the widths of Fe II 
and the VBLR (r = 0.66, P < 0.0001 for the Ha 
VBLR, and r=0.45, P < 0.0001 for H/3 VBLR). 

Relationships between the shifts of Fe II and 
Ha and H/3 components were also considered ( Ta- 
ble EJ. We found correlations with the shifts of the 
Ha ILR (r = 0.30, P = 0.0004) and the H/3 ILR 
component (r = 0.39, P < 0.0001), but no corre- 
lation between the shift of Fe II and the shifts of 
other Ha and H/3 components. 

We found that Fe II lines tend to have an av- 
eraged redshift of 270±180 kms -1 with respect to 
the transition wavelength, and with respect to nar- 
row lines 100±240 kms -1 (see Figure [16]) . 

3.3. Correlations between the Fe II line 
groups and other emission lines 

Correlations between the luminosity of the Fe II 
multiplet groups, total Fe II and luminosities of 
the [N II] and [O III] lines are presented in Table 
[51 In the L [N II] vs. L Fe II plot, seven outliers are 
observed with negligible [N II] . They are probably 
caused by an underestimate in the fit of the [N II] 
lines due to the blending with the Ha line. With- 
out these outliers, correlation is ~ 0.70, P<0.0001. 
There is no correlation between L [N II] /L [O III] 
vs. L Fe II. 

We analyzed relationships among the luminosi- 
ties of the Fe II line groups and the NLR, ILR 
and VBLR Ha components (Table®. All three 
Ha components are correlated with the Fe II line 



groups. Correlations are stronger with the ILR 
and VBLR components than with the NLR com- 
ponents. The same analysis was carried out for 
the luminosities of Fe II line groups and the NLR, 
ILR and VBLR components of H/3 (Table®. NLR 
component of H/3 also has a weaker correlation 
with Fe II lines (r~0.60, P<0.0001) than broad 
H/3 components (r~0.90, P<0.0001). 

We also investigated the anti-correlations be- 
tween EW Fe II and EW [O HI] and between EW 
Fe II and EW [O III]/EW H/3 (jBoroson fe Green 
1992). We confirmed the existence of these rela- 



tionships in our sample, with Fe II lines separated 
into F, S, and G l ine groups (Figure IT71 Ta ble l9l). 
A difference from iBoroson fe Greenl (|1992[ ) , who 
measured the equivalent widths of all lines refer- 
ring to the continuum level at the A4861 A, is that 
we calculated the equivalent widths by estimating 
the continuum below all the Fe II lines considered 
(see F igure El). Note also that IBoroson fe Green 
(|1992[ ) used mainly measurements from lines from 
multiplets 37 and 38. 

We found a correlation between total EW Fe II 
and EW [O III] r = -0.39, P < 0.0001. The 
correlation coefficient was significantly lower for 
the lines obtained from I Zw 1 (r = —0.20, P = 
0.0006) than for the other three groups (see Ta- 
ble [9]). We also analyzed relationships between 
EW Fe II line groups and EW [O III]/EW H/3. 
For total EW Fe II we found r= - 0.46, P<0.0001 
(Figure [T71 right), but for lines from I Zw 1 group 
r= - 0.28, P<0.0001. 

No significant correlations were found between 
EW Fe II and EW [N II] lines (Table E}. 

We also investigated whether there were any 
trends between equivalent widths of Fe II F, S, 
G line groups and equivalent widths of NLR, ILR 
and VBLR Ha and H/3 components (Table E}. We 
found no significant correlations. 

The correlations between EW Fe II and the 
Doppler widths of H/3 components are also con- 
sidered (Table [TO)) . As was expected, an inverse 
correlation was found between EW Fe II and the 
widths of the broad H/3 components (ILR and 

VBLR, as well as for the H/3 FWHM, ? 0.30, 

P < 0.0001). These correlations are part of EV 1. 
But, contrary to this, we found a positive correla- 
tion (r = 0.30, P < 0.0001) between the EW Fe II 
and the width of NLR component. 
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3.4. Relations among emission line strength 
and continuum luminosity 

In Sec 3.3. we have confirmed the anti- 
correlation EW Fe II vs. EW [O III] (also EW 
Fe II vs. EW [0 IH]/EW H/3). To try to under- 
stand the underlying physics which governs this 
anti-correlation, we examined its dependence of 
the continuum luminosity and redshift. 

Because of selection effects, cosmological red- 
shift, z, and continuum luminosity are highly cor- 
related in our sample, so it is difficult to clearly 
distinguish which influence dominates in some cor- 
relations. Since our sample has a more uniform 
redshift distribution than continuum luminosity 
distribution (see Figure [2J , we binned equivalent 
widths of the [O III] and Fe II lines within the 
z = - 0.7 range, using a Az = 0.1 bin size. The 
binned data are presented in Figure [18] It can be 
seen that as redshift increases in the z < 0.4 range, 
EW Fe II also increases, but EW [O III] decreases 
(Figure [18]) . For z > 0.4, the trend is not so obvi- 
ous, probably due to the larger scatter of the data 
at higher redshifts. 

We derived correlations coefficients for the ra- 
tios of the [O III], Fe II and H/3 lines and Lasioo 
(also with redshift, see Figure [T9l Table [TT]). A 
significant anti-correlation was found for the ra- 
tio of EW [O III]/EW Fe II vs. L A5 i 00 and for 
EW [O m]/EW Fe II vs. z (r = -0.46, P = 
and r = —0.48, P = 0, respectively). We found 
that the ratios of EW [O III]/EW and EW 
H/3 NLR/EW Fe II also decrease as La 5 i o (red- 
shift) increases. As in the previous case, the dif- 
ference between the P-values for correlation with 
La51oo and with z is very small, so we cannot tell 
which effect dominates. In contrast to this, for 
EW HAotal/EW Fe II vs. L A5 ioo, no significant 
trend is observed, but for EW H/3 tota i/EW Fe II 
vs. z, there is a weak correlation (r = 0.28, P = 
1.3E-6). 

Because of discrepancies in the properties ob- 
served between sub-samples within different H/3 
FWHM ranges (see Sec 3.1), we analyzed the re- 
lationships separately for the sub-samples with 
FWHM H^ <3000 kms^ 1 and FWHM H/3 >3000 
kms -1 . It will be noticed that, in general, all 
the correlations considered are stronger for the 
FWHM H/3 >3000 kms" 1 sub-sample than for 
spectra with narrower H/3 lines (see Table fTTj) ■ 



The only exception is the correlation of EW H/3 
NLR/EW Fe II vs. L>5ioo, which is more signifi- 
cant for the FWHM H/3 <3000 kms^ 1 subsample. 

It is interesting to connect these anti-correlations 
(EW Fe II vs. EW [O III] and consi dered ratios vs. 
LxRinn ) with the Baldwin effect. iBaskin fc Laoi 
( 20041 ) confirmed a correlation between the Bald- 



win effect and some of the emission parameters 
which define Eigenvector 1 correlations (which arc 
related to EW [O III] - EW Fe II anti-correlation). 
Also, it has been found that the [O III] lines show 
strong Baldwin effect ([Dietrich et al.1 120021 ). On 
the other hand, no Baldwin effect has been noticed 



for th e optical Fe II and H/3 lines ([Dietrich et al 



2002), or e ven an invers e Baldw in Effect has been 



reported: ICroom et al.l (120021) found an inverse 
Baldw in effect in H/3 and Netzer fc Trakhtenbrotl 
(|2007l ) found one for both H/3 and optical Fe II 
lines. 

Because of this, we examined the dependence 
of the equivalent widths of Fe II, H/3 and [O III] 
lines vs. L\^iqq and z in our sample. Since H/3 
is decomposed in NLR, ILR and VBLR and the 
iron lines are separated in the groups they can be 
considered in more detail. Also, we performed sep- 
arate analyses for the sub-samples with different 
H/3 FWHM ranges (Table [TT]). 

Analyzing the total sample (i.e., the whole H/3 
FWHM range), it is obvious that Fe II lines show 
an inverse Baldwin effect, which is specially strong 
for the central (r = 0.33, P < 0.0001) and red part 
(r = 0.44, P < 0.0001) of the Fe II shelf (S and G 
groups), as well as for the group of lines from I Zw 
1 (r = 0.32, P < 0.0001). The correlation of total 
Fe II vs. L A5 ioo is presented in Figure [201 What 
is interesting is that a strong inverse Baldwin ef- 
fect is not observed for the F lines (blue part) 
(r = 0.16, P = 0.004). In the relationship between 
H/3 components and continuum luminosity for to- 
tal sample, no trend is observed for broad H/3 (ILR 
and VBLR) , but it can be noticed that the narrow 
H/3 component anti-correlates with continuum lu- 
minosity (r = -0.36, P < 0.0001), - i.e. they 
show a Baldwin effect - as well as the [O III] lines, 
for which the previously found trend is confirmed 
(r = -0.43, P < 0.0001 - see Figured)])- 

If we compare these correlations with the cor- 
responding correlations with redshift, it is obvious 
that the iron lines (total Fe II and the multiplet 
groups) correlate more strongly with redshift than 
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with luminosity. In the case of other lines we con- 
sider, the difference between the correlations with 
Lasioo and z, is not so significant. Generally, it 
seems that the lines which show an inverse Bald- 
win effect (Fc II and H/3 VBLR) correlate more 
significantly with redshift than with luminosity, 
while the lines which show a (normal) Baldwin ef- 
fect ([O III] and H/3 NLR) have more significant 
correlations with continuum luminosity. 

We find that the coefficients of correlations be- 
tween the EWs of lines and the continuum lumi- 
nosity (and z) depend on the H/3 FWHM range 
of the sub-sample. It seems that continuum lumi- 
nosity affects iron lines and H/3 more for the sub- 
sample with narrower H/3 line (FWHM H/3 <3000 
kms -1 ). In that sub-sample all Fe II groups show a 
stronger inverse Baldwin effect, H/3 NLR decreases 
more strongly with an increase in continuum lu- 
minosity, and an inverse Baldwin effect is also ob- 
served for H/3 VBLR (r = 0.41, P < 0.0001). 
This is not seen in the FWHM H/3 >3000 kms" 1 
sub-sample. In contrast with this the Baldwin ef- 
fect is stronger in the [O III] lines for the FWHM 
H/3 >3000 kms -1 subsamplc. 

4. Discussion 

From our analysis of optical Fe II lines (AA 
4400-5400 A), we can try to investigate physical 
and kinematical characteristics of the Fe II emit- 
ting region in AGN as well as the connection be- 
tween the Fe II and other emission regions. We 
have investigated above correlations between the 
Fe II emission properties and some spectral fea- 
tures. Correlations of the Fe II lines were con- 
sidered separately for different multiplet groups, 
which enable more detailed analysis. Here we give 
some discussion of the results we have obtained. 

4.1. The Fe II line group ratios possible 
physical conditions in the Fe II emit- 
ting region 

Although we used very simple approximations 
to calculate the intensities of the observed Fe II 
lines, we found that the calculated intensities can 
satisfactorily fit the Fe II shelf within AA 4400- 
5500 A range. This approach is very simpli- 
fied, but in some cases, it enables better fit than 
much more complicated theoretical models (see 
Appendix B). We included excitation temperature 



in calculation and found that it is in the most 
cases around 10000 K (9646±2143 K), see Figure 
IH Roughly estimated temperature of the Fe II 
emitting region from our fitting (assumed tem- 
plate) is in a good agreement with t he predict ion 
earlier given in literature (-7000 K. |jolvlll987h . 

We found that the most intensive emission in 
the optical part arises in transitions with the lower 
term b 4 F (multiplets 37 and 38). We also con- 
sidered the ratios of multiplet groups. The line 
ratios may indicate some physical properties. For 
example, it has been shown that the Balmer line 
ratios are velocity dependent in AGN (see e.g . 



1986: lStiroe l fl990. 



I Shuder lll98lll984l:ICrenshaw 
ll99lHSnedden fe Gaskell 112007ft and this is proba 
bly related to a range of physical conditions (elec- 
tron temperature and dens ity) and to the radia - 
tive transfer effects (see e.g. lPopovic Il2003ll2006h . 
In general, the ratios between F, G and S groups 
indicate the ratio between the blue, red and cen- 
tral part of the Fe II shelf around H/3 lines. As 
can be seen in Table there is correlation be- 
tween F/G, F/S and G/S ratios vs. FWHM 
of H/3 line. It is interesting that this correlation 
is not present when we consider the cases where 
FWHM(H/3)< 3000 kms" 1 . We should note here 
that for the case FWHM(H/3)< 3000 kms -1 we 
had only AFWHM(H/3)«2000 km s" 1 and this 
may affect the obtained results. As we men- 
tioned above, the separation of the sample into 
two sub-samples a ccording to the H/3 F WHM is 
similar as given in ISulentic et al. (I2009L and ref- 
erence therein). It seems that the characteristics 
of the Pop A (FWHM H/3 <4000kms- 1 ) and B 
(FWHM H/3 >4000kms- 1 ) of AGN (as proposed 
by the above-mentioned authors) can be recog- 
nized in the Fe II line ratios: namely, the correla- 
tions and trends observed in Table [5] indicate that 
the ratios of blue and red (as well as blue and cen- 
tral) parts of Fe II shelf, anti-correlate with con- 
tinuum luminosity more strongly in Pop A than in 
Pop B. Also, for Pop B, all considered flux ratios 
(F/G, F/S and G/S) increase with increasing of 
H/3 width and with decreasing Fe II shift, while 
that trend is not observed in Pop A. 

The observed anti-correlations of F/G and F j S 
ratios vs. Lasioo could be affected by the differ- 
ing degrees of inverse Baldwin effect observed for 
the Fe II lines from the three line groups since the 
equivalent widths of the S and G groups increase 
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more significantly with L\ 5100 (r = 0.33 and 
r = 0.43, P < 0.0001, respectively), than those 
of the F group (r = 0.16, P = 0.04, see Table QT). 
Note here that the ratio of G/S does not correlate 
with the continuum. Also, since the inverse Bald- 
win effect for the S and G groups is stronger for 
the FWHM H/3 <3000 kms" 1 sub-sample, the cor- 
relations of F/G and F/S vs. Lasioo are stronger 
for FWHM H/3 <3000 kms" 1 one (Table ©. 

From the above-mentioned correlations, we can 
ask some intriguing questions about the physi- 
cal conditions and processes in the Fe II emit- 
ting region: (1) which processes cause the increase 
of Fe II equivalent width with increasing L>5ioo? 
Note that the equivalent widths of the majority 
of other emission lines in AGN spectra decreas e 
with increasing L A5100 (see iDietrich et al.l 120021) ; 
(2) why is the inverse Baldwin effect correlation 
lower for the F group comparing to the S and G 
groups? and (3) what causes that inverse Baldwin 
effect for Fe II (in S and G) to be more signif- 
icant in the sub-sample with FWHM H/3 <3000 
kms -1 , than in one with broader H/3 line? Answer- 
ing some of these questions may lead to a better 
understanding of the complex Fe II emitting re- 
gion. 

4.2. Location of the Fe II emitting region 



.Marziani fe Sulentic I (|1993[ ) and lPopovic et al.. 
( 20041 ) noted earlier that Fe II lines may originate 
in the intermediate-line region, which may be the 
transition from the torus to the BLR. Recently, 
other authors also found that the optical Fe II 
emission may arise from a region in the outer parts 
of the BLR or several times larger than the H/3 one 



( Hu et al.ll2008atlKuehn et al. 2008; [Popovic et al 
2004 l2007lT2009i: iGaskell et al.ll2007f ) 

To explore the kinematics of the Fe II emit- 
ting region, which may indicate the location of the 
Fe II emission, we compared the derived widths 
and shifts of the Fe II and the H/3 and Ha com- 
ponents. As can be seen in Figure [Ml and Figure 
H"5| there is an indication that Fe II emitting re- 
gion is located as the IL-emitting region. More- 
over, the correlation between the Doppler widths 
of the Fe II and IL regions is significantly higher 
than for the VBLR and NLR. There is no signif- 
icant correlation between the Doppler widths of 
the Fe II and NL regions (r = 0.01, P = 0.87 
for Ha, r = 0.05, P = 0.40 for H/3), but there 



is a correlation between the Fe II and VBLR re- 
gions (r = 0.66, P < 0.0001 for Ha, r = 0.45, 
P < 0.0001 for H/3). This may indicate that one 
more component of the Fe II lines is arising from 
the VBLR. Note here that, due to the complex 
Fe II template, we assumed only one component 
for each Fe II line (see Sect. 2.2). The reason- 
able fits of the Fe II shelf indicate that, the VBLR 
component of Fe II should be significantly weaker 
than the ILR one. 

These indications are also supported by the cor- 
relation between the shift of the Fe II lines and 
the ILR components of Ha and H/3. The distribu- 
tion of the Fe II line shifts is shown in Figure [16j 
The shift is obtained with respect to the transi- 
tion wavelength (Figure [16j left) and with respect 
to the narrow lines (Figure [16J right). We found 
that Fe II lines are slightly re dshifted, which i s 
in agreement with the results of Hu et al. ( 2008ah . 
But in this case, the average value of the Fe II shift 
relative to the transition wavelength is 270 ± 180 
kms -1 , while the average value of the Fe II shift 
with respect to the narrow lines is 100 ±240 kms" 1 . 
This im plies a significantl y smaller redshift than 
found in iHu et al.l (|2008al) . A slightly redshifted 
Fe II emission may indica te that Fe II l ines arise 
in an inflowing region fsee IGask ell 2009). 

4.3. EW Fe II vs. EW [O III] anti- 
correlation 

One of the problems mentioned in the introduc- 
tion is the anti-correlation between the equivalent 
widths of the [O III] and Fe II lines which is relate d 
to EV1 in the analysis of lBoroson fe Greenl (|1992l ) . 
Some physical causes proposed to explain Eigen- 
vector 1 correlations are: ( a) E ddington ratio 



L/L_Rdd (jBoroson fe Greenl Il992t iBaskin fe Laor 



2004), (b) black hole mass Mb h , and (c) inclina- 



tion a ngle dMilev fe Millerlll979l: IWills fe Browne! 



19861: iMarziani et al.l l200ll T~ fWang et al.l (|2006l) 

also suggested that EV1 may be related to AGN 
evolution. 

To try to understand the EW Fe II vs. EW 
[O III] anti-correlation, we examined its relation- 
ship to continuum luminosity and redshift (see 
Sec 3.4, Table fTTT) . We found an anti-correlation 
between the EW [O III]/EW Fe II ratio and 
La51oo- Also, we examined the relations of equiv- 
alent widths of Fe II and [O III] lines vs. L>5ioo. 
We confirmed a strong Baldwin effect for [O III] 
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lines, and an inverse Baldwin effect for EW Fe II 
lines, i.e. we found that as continuum luminosity 
increases, EW Fe II also increases, but EW [O III] 
decreases. This implies that the EW Fe II - EW 
[O III] anti-correlation may be influenced by Bald- 
win effect for [O III] and an inverse Baldwin effect 
for Fe II lines. Also, in our analysis we found 
that the strength of the Baldwin effect depends 
on the H/3 FWHM of the sample (see Table EJ. 
Note that H/3 FWHM is one of the parameters 
in Eigenvector 1. As it is mentioned in the in- 
troduction (Sec 1), some indications of connection 
be tween Baldwin e ffect a nd E V1 have been give n 



.Baskin fe Laorl (l2004h and bong etaD <|2009h . 
Ludwig et al.l (|2009l )~ found that the significance of 



the EV1 relationships is a strong function of con- 
tinuum luminosity, i.e., the relationships between 
EW Fe II, EW [O III] and H/3 FWHM can be de- 
tected only in a high-luminosity sample of AGNs 
(logAL 5 ioo > 44.7), as well as a Baldwin effect for 
[O III] linefS 

The origin of the Baldwin effect is still not un- 
derstood and is a matter of debate. The increase 
of the continuum luminosity may cause a decrease 
of the covering factor, or changes in the spec- 
tral energy distribution (softening of the ionizing 
continuum) which may result in the decrease of 
EWs. The inclination angle m ay also be related t o 
Baldwin effect (for review see iGreen et al . 2001). 
The physical properties which are usually consid- 
ered as a primary cause of the Baldwin ef fect are: 
M Rff d Warner et al .1120031: JXu et alj|2008h , L/LQ 
(|Boroson fe Greenl 1199a Baskin fe Laorl l2004h 
and changes in gas metalicity (jPietrich et al.l 
2002). Also, a connection between B aldwin effect 



and AGN evolution is possible (see IGreen et al 



2001) 



We investigated if the correlations observed be- 
tween Lasioo and the equivalent widths of the lines 
(as well as EW [O III] vs. EW Fe II) are primarily 
caused by evolution. Since continuum luminosity 
is strongly correlated with cosmological redshift 
in our sample, it is difficult to separate luminosity 
from evolutionary effect. The exception are corre- 
lations of EW Fe II vs. z, as well as the ratio EW 
H/?totai/EW Fe II vs. z, which are significantly 



stronger than correlations of the same quantities 
with Lasioo- This implies that inverse Baldwin ef- 
fect of Fe II may be governed first of all by an 
evolutionary e ffect. This r esult is in agreement 
with results of IGreen et al. ( 2001 ). 



7 The majority of objects in our sample have continuum 
luminosities in the range 44.5 <logAZ/5ioo < 46 range (see 
Figd. 



5. Conclusions 

In this paper we have investigated character- 
istic of the Fe II emission region, using the sam- 
ple of 302 AGN from SDSS. To analyze the Fe II 
emission, an Fe II template was constructed by 
grouping the strongest Fe II multiplets into three 
groups, according to the atomic properties of tran- 
sitions. 

We have investigated the correlations of these 
Fe II groups and their ratios with other lines in 
spectra. In this way, we tried to find some physical 
connection between the Fe II and other emitting 
regions, as well as to connect Fe II atomic struc- 
ture with the physical properties of the emitting 
plasma. 

We also investigated the kinematical connection 
between Fe II and the Balmer emission region. In 
particular the anti-correlation of EW Fe II - EW 
[O III] and its possible connection with AGN lu- 
minosity and evolution were analzed. From our 
investigation we can outline the following conclu- 
sions: 

1. We have proposed here an optical Fe II tem- 
plate for the AA4400-5500 A range, which 
consists of three groups of Fe II multiplets, 
grouped according to the lower terms of 
transitions (F, S and G), and an additional 
group of lines reconstructed from the I Zw 
1 spectrum. We found that template can 
satisfactorily fit the Fe II lines. In spectra 
in which Fe II emission lines have different 
relative intensities than in I Zw 1, this tem- 
plate fit better than empirical and theoreti- 
cal templates based on I Zw 1 spectrum (see 
more in Appendix B). Using this template, 
we are able to consider different groups of 
transitions which contribute to the blue, cen- 
tral and red part of Fe II shelf around the 
Rf3+[0 III] lines. 

2. We find that the ratios of different parts of 
the iron shelf (F/G, F/S, and G/S) depend 
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of some spectral properties such as: contin- 
uum luminosity, H/3 FWHM, shift of Fe II, 
and Ha /H/3 flux ratio. Also, it is noticed 
that spectra with H/3 FWHM greater and 
less than ^3000 kms -1 have different prop- 
erties which is reflected in significantly dif- 
ferent coefficients of correlation between the 
parameters. 

We found that the Fe II emission is mainly 
characterized with a random velocity of 
^1400 kms -1 that corresponds to the ILR 
origin, which is also supported by the sig- 
nificant correlation between the width of 
Fe II and Ha, H/3 ILR widths. This is 
in agr eement with the previous investiga- 
tions (IPopoyic et al.| [20o3 iHu et al.|[2008at 



Kuehn et al.l 120081 ). but unlike the earlier 
investigations, we found a slight correlation 
with the width of VBLR. Therefore, it is 
possible that Fe II is partly produced in the 
VBLR and cannot be resolved from contin- 
uum luminosity. Also, we found that the 
Fe II lines are slightly redshifted relative to 
the narrow lines (~ 100 kms -1 ) (see Figure 

The Balmer lines were decomposed into 
NLR, ILR and VBLR components and 
relationships among H/3 components and 
some spectral properties were investigated. 
We found a positive correlation between 
the width of narrow lines and EW Fe II 
(r = 0.30, P < 0.0001), while the width 
of the broad H/3 components anti-correlates 
with EW Fe II. Also, we found a Baldwin 
effect in the case of H/3 NLR component, 
while the ILR component showed no correla- 
tion with continuum luminosity, and VBLR 
component shows an inverse Baldwin effect 
for the FWHM < 3000 kms" 1 subsample. 

We confirm in our sample the anti-correlation 
between EW Fe II and EW [O I II] which is 
relate d to Eigenvector 1 (EV1) in lBoroson fe Greenl 
(|1992n and we examined its dependence on 
the continuum luminosity and redshift. We 
found an inverse Baldwin effect for Fe II 
lines from the central and red part of the 
Fe II shelf (S and G), but for Fe II lines 
from blue part (F) no correlation was ob- 
served. A Baldwin effect was confirmed for 



the [O III] lines. Since EW Fe II increases, 
and EW [O III] decreases with increases of 
continuum luminosity (and also the trend 
of decreasing of EW [O III]/EW Fe II ra- 
tio with luminosity is obvious) , the observed 
EW Fe II vs. EW [O III] anti-correlation 
is probably due to the same physical reason 
which causes the Baldwin effect. Moreover, 
it is observed that the coefficients of corre- 
lation due to Baldwin effect depend on H/3 
FWHM range of a sub-sample, which also 
implies the connection between the Bald- 
win effect and EV1. A decreasing trend is 
observed for EW [O III]/EW H/3 vs. contin- 
uum luminosity. 

6. We found that the equivalent width of the 
Fe II lines, as well as the EW H&otai/EW 
Fe II ratio have a stronger correlation with 
redshift, than with continuum luminosity. 
This implies that the inverse Baldwin effect 
of Fe II may be primarily caused by evolu- 
tion. 
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A. The sample selection 



Using an SQL search with requirements mentioned in Sec 2.1, we obtained 497 AGN spectra with broad 
emission lines. From that number, 188 spectra were rejected because of noise in the iron emission line region. 
An example of rejected spectrum is shown in Figure 1211 Also, three spectra were rejected because of bad 
pixels, and one because of very strong double-peak emission in H/3 (SDSS J154213. 91+183500.0). 

Three more spectra with broad H/3 and practically without the Fe II emission were rejected since it 
was not possible to fit the Fe II lines (SDSS J132756.15+111443.6, SDSS J104326.47+110524.2 and SDSS 
J083343. 47+074654. 5). These three spectra are shown in Figure |2"21 There is also a possibility that in 
these three objects the Fe II lines are very broad and weak so that they cannot be distinguished from the 
continuum emission. 

The sample contains 302 AGN spectra (table with the SDSS identification and obtained parameters from 
the best fit is available electronically only, as an additional file at http://www.aob.rs/~jkovacevic/table.txtl. 

To test the host galaxy contribution in the selected sample of 302 AGNs, we have comp ared our sample 
with the set of AGNs which host galaxy fraction is determined in lVanden Berk et al.l (]2006f ). We found 106 
common objects. In the case of 48 AGNs (45% of the sample) there is practically no host galaxy contribution 
(0%), and in the rest of the objects, the contribution of the host galaxy is mainly smaller than 20%. 



B. Comparison with other Fe II templates 

We compared two Fe II template s (one empi r ical a nd one theoretical) with our template. The empir- 
ical Fe II templat e was taken from iDong et al.l (|2008h using 46 broad and 95 narrow lines identified by 
Veron-Cettv et al.l ([20041 ) in the I Zw 1 spectrum within a AA4400-5500 A range. The Fe II lines were 
fit with 6 free parameters (the width, shift and intensity for narrow and the width shift and intensity 
for broad lines). We also considered the theoretical model from iBruhweiler fe Verner (|2008h calculated for 



log[nH/(cm~ 3 )]=ll, [£]/(l km s _1 )=20 and log[$H/(cm~ 2 s~ 1 )]=21, since we found that this gave the best 
fit for those values of physical parameters. With this model, Fe II lines were fit with 3 free parameters 
(width, shift and intensity). 

Figure [23] shows the spectrum of SDSS J020039. 15-084554.9, an object in which the iron lines have 
similar relative intensities to I Z w 1. We fit that spectrum with our template (a), with a template based 
on line intensities from I Z w 1 ( Dong et al. 20081 ) (b) and with a template calculated by CLOUDY code 
( Bruhweiler fe Verner 2008) (c). We found that for this object all three templates fit the iron lines well. 

In the case of the SDSS J141755.54 + 431155.8, however (Figure [24)) . the iron lines in the blue bump 
are slightly stronger compared to those in the red, and our model fit iron emission slightly better than the 
other two templates. The discrepancy between the blue and red Fe II bumps is especially emphasized in 
the spectrum of the object SDSS J111603. 13 + 020852.2 (Figure |25|). which makes a significant difference for 
comparison with the Fe II properties of I Zw 1 object. In this case, our template shows a small disagreement 
for the lines whose relative intensity is taken from I Zw 1, but other two models, cannot fit well this type of 
the Fe II emission. Disagreement is particularly strong in the blue bump (i.e., in the multiplets 37 and 38 of 
the F group). 



The relative intensities of Fe II lines 



To estimate the relative intensities of Fe II lines within the three grou ps, we stated tha t the intensit y 
of a line (m — > n) should be proportional to the number of emitters ( N m ) ( Oster brock 19891 : Griem 1997 ). 
i.e. for optically thin plasma (for Fe II lines) one can write: I mn = const ■ f N m dx, where A mn is the 
probability of the transition, X rnn is the transition wavelength, and £ is the depth of the emitting region. 

Of course, the density of emitters can be non-uniform across the region, but we can approximate here that 
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it is uniform, and in the case of non-thermodynamical equilibrium it could be written ( Osterbrock 19891) : 

N m ~ b(T, N e )g m exp (-E mn /kT), 

where b(T,iV e ) represents deviation from thermodyn amical equili brium. In the case of lines which have the 
same lower level, one may approximate the ratio as (|Griemlll997l) : 



h _ b i(T, N e ) ,A 2 , 3 /i 51 (Ei-E2)/kT 
h b2{T,N e y\i } h ' 92 ' 

Assuming that j^ e j w 1 we obtained Eq. (1) for estimation of line ratios within one group. 



(CI) 
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Fig. 1. — Grotrian diagram showing the strongest Fe II transitions in the AA 4400-5500 A region (top). Lines 
are separated into three groups according to the lower level of transition (middle) : F (dashed line) , S (solid 
line) and G (dash-dotted line). Bottom: the lines from the three line groups (solid line) and lines measured 
from I Zw 1, represented with dots. 
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Fig. 2. — Distribution of continuum luminosity (left) and cosmological redshift (right) in the AGN sample. 
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Fig. 3. — An example of removing the continuum from the spectrum of SDSS J075101.42 + 29174000.00; 
The level of the continuum is determined by interpolating between the chosen points at wavelengths without 
emission lines (top) . The spectrum without continuum emission is shown at the bottom; 
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Fig. 4. — The distribution of the excitation temperatures for the sample of 302 AGNs. Temperatures arc 
obtained by fits using Eq. (1). The mean value and dispersion are shown in the Figure. 
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Fig. 5.— Balmer lines in SDSS J111603. 13 + 020852.2 spectrum. Top: H/3 with a highly redshifted H/3 
VBLR component; bottom: Ha where the VBLR component is not redshifted. 
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Fe II (X 4924 A) 
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Wavelength (in A) 

Fig. 6.— The spectrum of SDSS J141755. 54 + 431155.8 in the AA 4400-5500 A region: (A) the observed 
spectra (dots) and the best fit (solid line). (B) H/3 fit with the sum of three Gaussians representing emission 
from the NLR, ILR and VBLR. The [O III] AA4959, 5007 A lines are fit with two Gaussians for each line of 
the doublet and He II A4686 A is fit with one broad Gaussian. The Fe II template is denoted with a dashed 
line, and also represented separately in panel C (bottom). 





350 r 




300 - 




250 - 




200 - 






c 

O 




c 

O 


150 - 


> 








o 


100 - 




50 - 




" 




-50 - 







: | 

Ha NLR 






[N II] 


\ | 
Ha ILR \ j 




Ha BLR \ I fl 







6500 



6550 6600 
Wavelength (in A) 



6650 



Fig. 7.— Fits of SDSS J141755.54+ 431155.8 in the AA 6400-6800 A region: Ha is fit with the sum of three 
Gaussians which represent emission from NLR, ILR and VBLR and [N II] AA6548, 6583 A lines are fit with 
one Gaussian for each line of doublet. Narrow [N II] lines are denoted with dashed lines. 
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Fig. 8. — Example of measuring of equivalent widths of the Fe II group of lines. Top: fit spectrum (SDSS 
J020039. 15-084554. 9); middle: Fe II lines and continuum; bottom: Fe II lines normalized on the continuum 
level. 
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Fig. 10. — Distribution of the ratios of luminosities of the three Fe II groups (F, S, and G). The mean 
value and the dispersion of the distribution are given in the plots. We indicate with a dashed vertical line 
the ratio of Fe II groups we obtained in the I Zw 1 spectrum. 
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Fig. 11 — The distribution of the ratios of total Fe II (4400-5500 
mean value and the dispersion of the distribution are shown. 
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Fig. 12.— The relationship between F/S and H/3 FWHM, for total sample (left), sub-sample of H/3 FWHM 
< 3000 kms- 1 (middle) and the sub-sample of R/3 FWHM > 3000 kms" 1 (right). 
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Fig. 13. — The relationship between Fe II group ratios and the Balmer decrement. 
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Fig. 14. — The widths of the Fe II lines compared to the widths of the H/3 (left) and Ha (right) components. 
On the X-axis are the widths of Fe II, and on Y-axis are the widths of the NLR (triangles), ILR (circles), and 
VBLR (squares) components of H/3 (Ha). The dotted vertical line shows the average value of Fe II widths, 
while the dotted horizontal lines show the average values of the H/3 (Ha) components. The average value of 
for the Fe II lines is the same as or very close to the average width value of the ILR components of H/3 and 
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Fig. 15. — The correlation between the widths of Fe II and the H/3 ILR (left), Ha ILR (middle) and 
Ha VBLR components (right). In all cases, correlations are observed (r = 0.67, 0.72, 0.62), indicating 
kinematical connection between these emission regions. 
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Fig. 16. — Distribution of the Fe II shift with respect to the transition wavelength (left) and distribution of 
the Fe II shift with respect to the [O III] lines (right). 
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Fig. 17.— Relationship between the EW [O III] A5007 A vs. EW Fe II (left) and EW [O III]/EW H/3 vs. 
EW Fell (right). 
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Fig. 18. — Relationship between the equivalent widths, EW [O III] A 5007 A and EW Fe II, binned for 
cosmological redshift. The error bars are the dispersions in each sub-sample within different redshift bins. 
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Fig. 19. — First panel: EW [O III]/EW Fe II ratio vs. continuum luminosity (AL5100); second panel: the 
same for EW [O III]/EW H/3 ratio. Objects with redshift within range z<0.1 are denoted with open squares, 
0.1<z<0.2 with filled triangles, 0.2<z<0.3 with open circles, 0.3<z<0.4 with filled squares, 0.4<z<0.5 with 
open triangles, 0.5<z<0.6 with filled circles and 0.6<z<0.7 with stars. [O III] includes only the A5007 A 
component. 
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Fig. 20. — The Baldwin effect is significant for the [O III] lines (first panel), while an inverse Baldwin effect 
is detected for the optical Fe II lines (second panel). Symbols have the same meaning as in Figure HU 
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Fig. 21. — Example of a spectrum rejected because of strong noise which inhibited a precise fit to the iron 
lines. 
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Fig. 22. — The three spectra with broad H/3 line and without Fe II emission. 
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Fi g. 23. — Examples of fits to SDSS J020039. 15-084554.9: with o ur template (a), with t he em pirical template 
of Dong et al. ( 20081 ) (b), and with the theoretical template of iBruhweiler fc Verner ( 20081 ) (c). Since this 
AGN has iron emission equally strong in blue and red iron bump (as I Zw 1), all three models fit the observed 
lines well. 
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Fig. 24. — Examples o f fits t o SDSS J141755. 54 + 431155.8: with our template (a), with th e emp irical 
template of Dong et al. ( 20081 ) (b), and with the theoretical template of Bruhweiler fc Verner ( 20081 ) (c). 
Iron emission is a bit stronger in the blue bump, than in red. Our model fits this spectrum slightly better 
than other two. 
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Fig. 25.— Examples o f fits t o SDSS J111603.13 + 020852.2: with o ur template (a), with the e mpirical 
template of iDong et al. ( 20081) (b) and with the theoretical template of Bruhweiler fe Verner ( 2008 ) (c). In 
this object iron emission is much stronger in blue than in red bump. Our template show disagreement for 
lines whose relative intensity is taken from I Zw 1, but other the three wavelength regions based on our three 
line groups fit the observed Fe II well. The other two models cannot fit this kind of Fe II emission well. 
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Table 1: A list of the 35 strongest Fe II lines within the A A 4400-5500 A region used to calculate the Fe II 

temp late. In the first column are wavelengths (in air), in the second multiplet numbers (jGiridhar fe Ferrol 

1995), in the third are terms of transitions, in the fourth are gf values used for the template calculation, 

in the fifth are the references for the source of oscillator strengths, and in the 6th-8th columns are relative 
intensities, calculated for T=5000 K, 10000 K and 15000 K. Intensities of lines from the F, S and G groups 

are normalized to intensities of the A4549.474 A, A5018.44 A, and A5316.6 A lines (respectively). 

Wavelength Multiplet Transition gf ref. Relative intensity 

T=5000 K T=10000 K T=15000 K 

4472.929 37 b 4 F 5/2 - z 4 F° /2 4.02E-04 1 O033 O036 O037 

4489.183 37 b 4 F 7/2 - z 4 F° /2 1.20E-03 1 0.105 0.110 0.111 

4491.405 37 b 4 F 3/2 - z i F° /2 2.76E-03 1 0.226 0.243 0.249 

4508.288 38 b 4 F 3/2 - z 4 D° /2 4.16E-03 2 0.345 0.367 0.375 

4515.339 37 b 4 F 5/2 - z 4 F° /2 3.89E-03 2 0.333 0.348 0.353 
4520.224 37 b 4 F 9/2 - z 4 F° /2 2.50E-03 3 0.235 0.233 0.233 
4522.634 38 b 4 F 5/2 - z 4 D° /2 9.60E-03 3 0.827 0.859 0.870 
4534.168 37 b 4 F 3/2 - z 4 F° /2 3.32E-04 3 0.028 0.029 0.030 
4541.524 38 b 4 F 3/2 - z 4 D° /2 8.80E-04 3 0.075 0.078 0.079 
4549.474 38 b 4 F 7/2 - z 4 D° /2 1.10E-02 4 1.000 1.000 1.000 
4555.893 37 b 4 F 7/2 - z 4 F° /2 5.20E-03 3 0.477 0.474 0.474 

4576.340 38 b 4 F 5/2 - z 4 D° /2 1.51E-03 4 0.136 0.136 0.136 
4582.835 37 b 4 F 5/2 - z 4 F° /2 7.80E-04 3 0.070 0.070 0.070 
4583.837 38 b 4 F 9/2 - z 4 D° /2 1.44E-02 2 1.420 1.353 1.331 
4620.521 38 b 4 F 7/2 - z 4 D° /2 8.32E-04 4 0.080 0.076 0.075 
4629.339 37 b 4 F 9/2 - z 4 F° /2 4.90E-03 4 0.497 0.459 0.447 
4666.758 37 b 4 F 7/2 - z 4 F° /2 6.02E-04 4 0.060 0.055 0.054 
4993.358 36 b 4 F 9/2 - z 6 P° /2 3.26E-04 4 0.041 0.030 0.027 

5146.127 35 b 4 F 7/2 - z 6 F° /2 1.22E-04 5 0.016 0.011 0.010 

4731.453 43 a °S 5/2 - z 4 fl° /2 1.20E-03 2 0.025 0.030 0.032 

4923.927 42 a 6 S 5/2 - z 6 P° /2 2.75E-02 4 0.656 0.693 0.706 

5018.440 42 a 6 S 5/2 - z 6 P° /2 3.98E-02 4 1.000 1.000 1.000 

5169.033 42 a 6 S 5/2 - z 6 P° /2 3.42E-02 1 0.929 0.854 0.831 

5284.109 41 a 6 S 5/2 - z 6 F° /2 7.56E-04 2 0.022 0.019 0.018 

5197.577 49 a 4 G 5/2 - z 4 F 3 ^ 2 7.92E-03 5 0.532 0.620 0.652 

5234.625 49 a 4 G 7/2 - z 4 F° /2 8.80E-03 3 0.615 0.695 0.723 

5264.812 48 a 4 G 5/2 - z 4 D° /2 1.08E-03 1 0.075 0.084 0.087 

5276.002 49 a 4 G 9/2 - z 4 F° /2 1.148e-02 2 0.861 0.928 0.951 

5316.615 49 a 4 G 11/2 - z 4 F° /2 1.17E-02 4 1.000 1.000 1.000 

5316.784 48 a 4 G 7/2 - z 4 D° /2 1.23E-03 5 0.089 0.097 0.099 

5325.553 49 a 4 G T/2 - z 4 F° /2 6.02E-04 4 0.044 0.047 0.048 

5337.732 48 a 4 G 5/2 - z 4 D° /2 1.28E-04 5 0.009 0.010 0.010 

5362.869 48 a 4 G 9/2 - z 4 D° /2 1.82E-03 5 0.142 0.146 0.148 

5414.073 48 a 4 G 7/2 - z 4 D° /2 1.60E-04 5 0.012 0.012 0.013 

5425.257 49 a 4 G 9/2 - z 4 F° /2 4.36E-04 5 0.035 0.035 0.035 



References. — (1) iFuhr et al.l (119811) , (2) iGiridhar fc Ferrol (119951 1. (3) NIST Atomic Spectra Database 
jhttp://physics.nist.gov/Ph ysRefData/ASD/ ), (4) iKuruczl ||1990D and (5) http:// www. pmp . uni-hannover . de/cgi- 
bin/ssi/tcst /kurucz / sckur . html 
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Table 2: List of the lines taken from Kurucz database ( [http://www.prnp.uni-hannover.de/cgi-bin/ssi/test/kurucz/sekur.html l. 
In first column are wavelengths (in air), in the second oscillator strengths, and in the third relative intensities 
measured in I Zw 1. 



Wavelength Transitions gf Relative intensity 



4418.957 y 4 G 5/2 - e 4 F 3/2 1.45E-02 31)0 

4449.616 y 4 G 9/2 - e 4 F 7/2 2.58E-02 1.50 

4471.273 y 4 D 5/2 - f 4 D 3/2 6.40E-03 1.20 

4493.529 y 4 G 10/2 - e 4 F 9/2 3.74E-02 1.60 

4614.551 y 6 P 5/2 - 6 D 7/2 2.69E-03 0.70 

4625.481 x 4 D 1/2 - f 4 D 1/2 8.51E-03 0.70 

4628.786 x 4 D 5/2 - f 4 D 5/2 1.83E-02 1.20 

4631.873 x 4 D 3/2 - f 4 D 3/2 1.34E-02 0.60 

4660.593 d 5 4s 2 2 G 7/2 - (4G)sp 2 G 7/2 1.15E-03 1.00 

4668.923 x 4 D 5/2 - f 4 D 7/2 3.13E-03 0.90 

4740.828 y 6 P 5/2 - e 6 F 7/2 1.98E-03 0.50 

5131.210 y 4 P 5/2 - e 4 D 5/2 2.13E-03 1.1 

5369.190 e 4 D 5/2 - 4 D 5/2 1.23E-03 1.45 

5396.232 x 4 F 5/2 - f 4 D 3/2 2.80E-03 0.40 

5427.826 b 4 G 11/2 - w 4 F 9/2 2.17E-02 1.40 



Table 3: Correlations between widths, shifts and luminosities of the NLR, ILR and VBLR components of the 
Ha and H/3 lines. The values in brackets are for the sample without the 16 objects with highly redshifted H/3 
VBLR components (>4000 kms -1 ). The coefficient of correlation (r) is in the second column, the P-value 
which is a measure of significance of the correlation is given in the third column, and in the next columns 
are the coefficients A and B from the fit with a linear function: H/3 = A ± B ■ Ha. The slopes are from 
regressions where X is taken as the independent variable. 



correlation (Ha and H/3) 


r 


P 


A 


B 


Ha NLR width vs. R0 NLR width 
Ha ILR width vs. H/3 ILR width 
Ha VBLR width vs. H/3 VBLR width 


0.99 
0.91 
0.50 (0.42) 


< 0.0001 

< 0.0001 

< 0.0001 


6.19 ±2.81 
165.37 ± 43.29 
1863.15 ±352.50 


0.97 ±0.01 
0.74 ±0.03 
0.54 ±0.08 


Ha NLR shift vs. H/3 NLR shift 
Ha ILR shift vs. H/3 ILR shift 
Ha VBLR shift vs. H/3 VBLR shift 


0.97 
0.90 
0.18 (0.41) 


< 0.0001 

< 0.0001 
0.041 (< 0.0001) 


16.40 ± 3.58 
-1.86 ± 16.97 
-93.64 ± 45.30 


0.91 ± 0.02 
1.09 ±0.05 
0.06 ±0.03 


L Ha NLR vs. L H/3 NLR 

L Ha ILR vs. L H(3 ILR 

L Ha VBLR vs. L H^ VBLR 


0.93 
0.96 
0.96 


< 0.0001 

< 0.0001 

< 0.0001 


3.40 ± 1.29 
1.93 ± 1.08 
0.70 ± 1.03 


0.93 ± 0.03 
0.97 ±0.03 
0.99 ±0.02 



Table 4: Correlatios between equivalent widths (EW) of the Fe II F, S and G multiplet groups. Relationships 
are fitted with function Y = A ± B ■ X . The coefficient of correlations, r, the corresponding measure of the 
significance of correlations, P, as well as the A and B coefficients, are shown in the table. 



X 


Y 


r 


P 


A 


B 


EW Fe II F 
EW Fe II F 
EW Fe II S 


EW Fe II S 
EW Fe II G 
EW Fe II G 


0.76 
0.80 
0.83 


< 0.0001 

< 0.0001 

< 0.0001 


0.21 ± 1.00 
-0.17 ± 1.29 
5.39 ±0.95 


0.66 ±0.03 
0.96 ± 0.04 
0.14 ±0.04 
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Table 5: Correlations between the ratio of Fe II line groups and spectral properties. In the first column is 
FWHM of broad H/3 (ILR+VBLR), in the second and third columns are the Doppler widths of ILR and 
VBLR H/3 components, in the fourth and fifth the Doppler width and shift of Fe II lines, and in the eighth 
column is the continuum luminosity at A5100 A. Correlation coefficients are calculated for total sample, for 
objects with H/3 FWHM < 3000 kms" 1 (1), and with H/3 FWHM > 3000 kms^ 1 (2). The correlations with 
P < 0.0001 are printed in bold print. 







H/3 FWHM 


w H/3 ILR 


w H/3 VBLR 


w Fe II 


sh Fe II 


log AZ/5100 


F/G 


r 


0.28 


0.24 


0.29 


0.18 


- 0.27 


- 0.32 




P 


7.13E-7 


2.49E-5 


2.37E-7 


0.002 


2.75E-6 


1.82E-8 


(1) F/G 


r 


- 0.10 


- 0.02 


0.16 


0.001 


- 0.10 


- 0.51 




P 


0.20 


0.831 


0.046 


0.987 


0.189 


5.72E-12 


(2) F/G 


r 


0.36 


0.27 


0.31 


0.17 


- 0.30 


- 0.36 




P 


1.22E-5 


0.001 


1.82E-4 


0.043 


3.31E-4 


8.98E-6 


F/S 


r 


0.41 


0.26 


0.16 


0.20 


- 0.27 


- 0.15 




P 


1.95E-13 


4.05E-6 


0.004 


5.47E-4 


2.44E-6 


0.011 


(1) F/S 


r 


- 0.21 


- 0.18 


0.07 


-0.23 


- 0.11 


- 0.41 




P 


0.009 


0.027 


0.407 


0.004 


0.182 


7.90E-8 


(2) F/S 


r 


0.59 


0.33 


0.12 


0.25 


- 0.29 


- 0.18 




P 


1.35E-14 


5.08E-5 


0.157 


0.002 


5.29E-4 


0.035 


G/S 


r 


0.26 


0.16 


0.03 


0.07 


- 0.11 


0.06 




P 


3.92E-6 


0.004 


0.646 


0.213 


0.048 


0.312 


(1) G/S 


r 


- 0.09 


- 0.17 


- 0.09 


- 0.26 


0.003 


0.14 




P 


0.243 


0.036 


0.277 


0.001 


0.965 


0.076 


(2) G/S 


r 


0.44 


0.28 


0.01 


0.15 


- 0.13 


0.006 




P 


6.15E-8 


8.26E-4 


0.867 


0.07 


0.113 


0.935 



Table 6: The relationship between log (L Ha/L H/3) and the ratios of the Fe II groups. 





log (L Ha/L H/3) 


log (F/G) 


r 


-0.36 




P 


< 0.0001 


log (F/S) 


r 


-0.34 




P 


< 0.0001 


log (G/S) 


r 


0.004 




P 


0.962 



37 



Table 7: Relationships between the Fe II width and the widths of the NLR, ILR, and VBLR components of the 
Ha and H/3 (left), and the same but for the shifts (right). Relationships are fitted with function Y = A+B-X. 
The coefficient of correlations, r, the corresponding measure of the significance of correlations, P, as well as 
the A and B coefficients, are shown in the table. 







w Fc II 


Corr. 






shift Fc II 


Corr. 


w Ha NLR 


r 


0.01 


No 


sh Ha NLR 


r 


0.02 


No 




P 


0.871 






P 


0.84 






A 


232.71 ± 27.12 






A 


149.78 ± 8.88 






B 


0.00 ± 0.02 






B 


0.00 ± 0.02 




w Ha ILR 


r 


0.77 


Yes 


sh Ha ILR 


r 


0.30 


Yes, weak 




P 


< 0.0001 






P 


0.0004 






A 


-94.35 ± 111.64 






A 


66.10 ± 26.11 






B 


0.96 ± 0.08 






B 


0.26 ± 0.07 




w Ha VBLR 


r 


0.66 


Yes 


sh Ha VBLR 


r 


0.01 


No 




P 


< 0.0001 






P 


0.930 






A 


609.84 ± 356.76 






A 


-49.14 ± 44.47 






B 


2.63 ± 0.26 






B 


0.01 ± 0.12 




w H/3 NLR 


r 


0.05 


No 


sh H,3 NLR 


r 


0.00 


No 




P 


0.401 






P 


0.928 






A 


275.84 ± 30.17 






A 


165.22 ± 9.63 






B 


0.02 ± 0.02 






B 


0.00 ± 0.03 




w H/3 ILR 


r 


0.73 


Yes 


sh H/3 ILR 


r 


0.39 


Yes 




P 


< 0.0001 






P 


< 0.0001 






A 


-114.23 ± 94.28 






A 


55.32 ± 14.80 






B 


1.17 ± 0.06 






B 


0.32 ± 0.04 




w H/3 VBLR 


r 


0.45 


Yes 


sh H/3 VBLR 


r 


-0.17 


No 




P 


< 0.0001 






P 


0.003 






A 


2241.03 ± 251.81 






A 


1301.87 ± 96.29 






B 


1.48 ± 0.17 






B 


-0.83 ± 0.28 





38 



Table 8: Correlations between luminosities (L) of the Fe II lines (total Fe II, F, S and G line groups) and 
luminosities of [N II], [O III] and Ha and H/3 components. Relationships are fitted with function Y = A+B-X. 
The coefficient of correlations, r, and the corresponding measure of the significance of correlations, P are 
shown in the table. 







log(L Fe II total) 


log(L F) 


log(L S) 


log(L G) 


log(L I Zw 1 group) 


log(L [N II]) 


r 


0.73 


0.72 


0.67 


0.73 


0.72 




P 


< 0.0001 


< 0.0001 


< 0.0001 


< 0.0001 


< 0.0001 


log(L [N II]/L [0 III]) 


r 


0.06 


0.08 


0.08 


0.10 


0.02 




P 


0.51 


0.38 


0.35 


0.25 


0.84 


log(L Ha NLR) 


r 


0.76 


0.74 


0.71 


0.75 


0.77 




P 


< 0.0001 


< 0.0001 


< 0.0001 


< 0.0001 


< 0.0001 


log(L Hq ILR) 


r 


0.91 


0.90 


0.84 


0.89 


0.90 




P 


< 0.0001 


< 0.0001 


< 0.0001 


< 0.0001 


< 0.0001 


log(L Ha VBLR) 


r 


0.90 


0.89 


0.82 


0.87 


0.90 




P 


< 0.0001 


< 0.0001 


< 0.0001 


< 0.0001 


< 0.0001 


log(L H/3 NLR) 


r 


0.64 


0.63 


0.61 


0.63 


0.64 




P 


< 0.0001 


< 0.0001 


< 0.0001 


< 0.0001 


< 0.0001 


log(L H/3 ILR) 


r 


0.93 


0.93 


0.88 


0.92 


0.92 




P 


< 0.0001 


< 0.0001 


< 0.0001 


< 0.0001 


< 0.0001 


log(L H/3 VBLR) 


r 


0.95 


0.95 


0.90 


0.93 


0.95 




P 


< 0.0001 


< 0.0001 


< 0.0001 


< 0.0001 


< 0.0001 



Table 9: The same as in Table 8, but for EWs. 







EW Fe II total 


EW F 


EW S 


EW G 


EW I Zw 1 group 


EW [N II] 


r 


0.03 


-0.00 


0.05 


0.07 


-0.01 




P 


0.74 


0.98 


0.58 


0.40 


0.89 


EW [O III] 


r 


-0.39 


-0.40 


-0.37 


-0.42 


-0.20 




P 


< 0.0001 


< 0.0001 


< 0.0001 


< 0.0001 


0.0006 


EW [O III]/EW H/3 


r 


-0.46 


-0.47 


-0.43 


-0.45 


-0.28 




P 


< 0.0001 


< 0.0001 


< 0.0001 


< 0.0001 


< 0.0001 


EW Ha NLR 


r 


-0.02 


-0.06 


-0.01 


0.01 


0.00 




P 


0.84 


0.49 


0.90 


0.89 


0.99 


EW Ha ILR 


r 


-0.15 


-0.12 


-0.12 


-0.18 


-0.13 




P 


0.08 


0.15 


0.15 


0.04 


0.14 


EW Ha VBLR 


r 


-0.17 


-0.19 


-0.18 


-0.19 


-0.03 




P 


0.05 


0.03 


0.03 


0.03 


0.69 


EW H/3 NLR 


r 


0.01 


0.01 


-0.09 


-0.04 


0.15 




P 


0.90 


0.91 


0.11 


0.46 


0.01 


EW H/3 ILR 


r 


-0.02 


0.02 


0.07 


-0.07 


-0.08 




P 


0.70 


0.69 


0.25 


0.24 


0.18 


EW H/3 VBLR 


r 


0.12 


0.11 


0.11 


0.04 


0.21 




P 


0.03 


0.07 


0.07 


0.52 


0.0003 
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Table 10: The equivalent width of Fe II vs. the widths of the H/3 components. 





EW Fe II 


width H/3 NLR 


r 0.30 
P < 0.0001 


width H/3 ILR 


r -0.31 
P < 0.0001 


width H/3 VBLR 


r -0.34 
P < 0.0001 


H/3 FWHM 


r -0.30 
P < 0.0001 



Table 11 

Correlations for the total sample (first and second column), the sub-sample with FWHM 
Hj3 <3000 kms" 1 (1), and sub-sample with FWHM H(3 >3000 kms -1 (2). The correlations 

WITH P<0.0001 ARE GIVEN IN BOLD PRINT. 







log AL5100 


log z 


(1) log AL 510 o 


(1) log z 


(2) log AL5100 


(2) log z 


log (EW [O III]/EW Fe II) 


r 


-0.46 


-0.48 


-0.44 


-0.45 


-0.55 


-0.56 




P 








5.67E-9 


1.89E-9 


1.69E-12 


6.40E-13 


log (EW [O III]/EW H/3) 


r 


-0.49 


-0.48 


-0.43 


-0.43 


-0.50 


-0.49 




P 








1.66E-8 


2.00E-8 


2.14E-10 


5.85E-10 


log (EW H/3 tota i/EW Fc II) 


r 


-0.19 


-0.28 


-0.28 


-0.31 


-0.45 


-0.49 




P 


0.001 


1.30E-6 


4.00E-4 


5.47E-5 


2.40E-8 


9.53E-10 


log (EW H/3 NLR/EW Fc II) 


r 


-0.42 


-0.42 


-0.52 


-0.51 


-0.31 


-0.31 




P 


2.35E-14 


4.75E-14 


2.42E-12 


7.17E-12 


2.22E-4 


1.97E-4 


log (EW [O III]) 


r 


-0.43 


-0.42 


-0.33 


-0.32 


-0.54 


-0.53 




P 


4.00E-15 


1.60E-14 


1.89E-5 


4.16E-5 


3.35E-12 


1.37E-11 


log (EW FeII tota i) 


r 


0.30 


0.39 


0.47 


0.53 


0.36 


0.41 




P 


1.89E-7 


3.82E-12 


7.05E-10 


6.09E-13 


1.06E-5 


4.16E-7 


log (EW Fe II F) 


r 


0.16 


0.26 


0.23 


0.28 


0.27 


0.35 




P 


0.004 


5.99E-6 


0.003 


2.81E-4 


0.001 


2.36E-5 


log (EW Fe II S) 


r 


0.33 


0.40 


0.52 


0.58 


0.37 


0.39 




P 


3.87E-9 


1.04E-12 


1.87E-12 


8.88E-16 


4.57E-6 


1.37E-6 


log (EW Fc II G) 


r 


0.44 


0.50 


0.61 


0.65 


0.50 


0.51 




P 


1.78E-15 











2.52E-10 


9.89E-11 


log (EW Fc II I Zw 1 group) 


r 


0.32 


0.23 


0.31 


0.39 


0.13 


0.18 




P 


1.77E-8 


4.58E-5 


8.33E-5 


4.29E-7 


0.120 


0.031 


log (EW H/3 NLR) 


r 


-0.36 


-0.33 


-0.44 


-0.41 


-0.21 


-0.20 




P 


2.05E-10 


6.78E-9 


8.87E-9 


8.79E-8 


0.010 


0.015 


log (EW H/3 ILR) 


r 


0.08 


0.03 


-0.04 


-0.04 


-0.11 


-0.08 




P 


0.161 


0.545 


0.611 


0.603 


0.21 


0.315 


log (EW H/3 VBLR) 


r 


0.16 


0.17 


0.41 


0.44 


-0.09 


-0.09 




P 


0.006 


0.002 


6.98E-8 


9.53E-9 


0.292 


0.271 
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